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We point out that the chiral Lagrangian describing pseudo-Goldstone bosons in partially quenched 
QCD has one more four-derivative operator than that for unquenched QCD with three flavors. The 
new operator can be chosen to vanish in the unquenched sector of the partially quenched theory. 
Its contributions begin at next-to-leading order in the chiral expansion. At this order it contributes 
only to unphysical scattering processes, and we work out some examples. Its contributions to 
pseudo-Goldstone properties begin at next-to-next-to-leading order, and we determine their form. 
We also determine all the zero and two derivative operators in the 0(p 6 ) partially quenched chiral 
Lagrangian, finding three more than in unquenched QCD, and use these to give the general form 
of the analytic next-to-next-to-leading order contributions to the pseudo-Goldstone mass and decay 
constant. We discuss the general implications of such additional operators for the utility of partially 
quenched simulations. 



I. INTRODUCTION 

Chiral perturbation theory (xPT) allows analytic calculations of low-energy QCD processes, results of which are 
given in terms of a number of undetermined constants, e.g. the Gasser-Leutwyler coefficients. These low-energy 
constants are fundamental QCD parameters that govern physical properties such as masses and scattering amplitudes. 
Lattice QCD provides a method for determining them from first principles, as long as simulations are done at small 
enough quark masses that X-PT (typically at next-to-leading order) is a good approximation. In practice, however, 
simulating light dynamical quarks is computationally expensive, and there is a limited range of quark masses where 
both lattice simulations are feasible and xPT is applicable. 

This situation can be improved using the partially quenched (PQ) approximation, in which valence quarks (those 
which appear in external states) and sea quarks (those which appear in loops) are allowed to have different masses. 
To ensure that valence quarks exist only in external states, valence quarks have bosonic ghost quark partners of 
equal mass. The presence of ghosts means that the PQ theory as a whole is unphysical, although the sea quark 
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FIG. 1: Parameter space of partially quenched theories. Physical theories, such as QCD, live on the diagonal line. Lattice 
simulations have been done with "light" quark masses as low as ~ m s trange/8. Note that this plot is schematic, since there are 
actually multiple valence and sea quarks, all of which can have different masses. 
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sector contained within it is physical. Consideration of PQ theories expands the parameter space available for lattice 
simulations and comparison to xPT, as shown in Fig. ^ 

In Refs. 0, 0] it is argued that PQ QCD can be used for a quantitative determination of the QCD low-energy 
constants. The primary assumption is that a generalization of xPT can be used to describe the low-energy behavior 
of PQ QCD. This theory, called PQxPT, must contain all of the operators of unqucnched xPT in order to describe the 
behavior of the unquenched sector. Quark mass dependence in both x?T and PQxPT is explicit, so the coefficients 
in the Lagrangian depend only on the number of quarks. Therefore the coefficients of operators in the unquenched 
chiral Lagrangian are identical to those of corresponding operators in the PQ chiral Lagrangian for N sea = 3. 1 

To make use of this observation one fits results from PQ simulations to the predictions of PQxPT and thereby 
determines the low-energy constants. Particular quantities studied in Refs. 0, y| were next-to-leading order (NLO) 
PGB masses and decay constants, as well as the hairpin vertex. It was asserted for these calculations that there 
are no new operators through NLO in the PQ chiral Lagrangian, so that no additional low-energy constants are 
needed at this order. Here we show that Refs. missed one operator in the NLO Lagrangian, and explore 
the consequences. While the existence of additional operators does not fundamentally change the approach, it does 
necessitate an elaboration of the program of using PQ QCD to determine low energy constants. 

The general situation is as follows. The PQ chiral Lagrangian contains two types of operators: first, those which, 
in the sea quark sector, reduce to operators present in the QCD chiral Lagrangian, and whose coefficients are thus 
identical to those in QCD, and, second, operators which vanish in the sea quark sector. We call the latter unphysical 
operators. Their coefficients are additional unphysical constants needed to describe the low-energy behavior of the 
PQ theory. They appear first at NLO in chiral perturbation theory, with increasing numbers required as one works 
at higher order. In general, both types of operator need to be taken into account when matching forms calculated in 
PQxPT to PQ lattice data. It turns out, however, that for masses, decay constants and the hairpin vertex, the new 
operators do not contribute until next-to-next-to-leading order (NNLO). Thus the results of Refs. [HQ are unchanged. 

This paper constitutes a preliminary investigation of the impact of the unphysical operators on the utility of PQ 
QCD. It is organized as follows. In Sec. [HI we derive the existence and form of the new, linearly-independent four- 
derivative operator present in the PQ chiral Lagrangian at NLO. We explore the consequences of this operator in 
Sec. Illll bv discussing the type of meson scattering processes to which it contributes. In Sec. lIVI we continue the study 
of Refs. 0, by calculating the leading contribution of this operator to meson masses and decay constants, which 
appear first at NNLO. This raises the general issue of extending the approach of using PQxPT to NNLO. In Sec.lVlwe 
take another step in this direction by discussing the analytic NNLO corrections to meson masses and decay constants. 
We determine the operators that contribute, construct a minimal, linearly-independent set, and show which of them 
vanish in the unquenched sector. Using these operators, we determine the general form of NNLO analytic meson mass 
and decay constant corrections in PQxPT, a result that is needed in fits to present lattice data. Finally, in Sec. IVII 
we summarize our findings, and discuss their general implications concerning the utility of PQ simulations. 

II. NEW FOUR-DERIVATIVE OPERATOR IN PQxPT 

We first recall the chiral effective Lagrangian describing the properties of the light pseudo-Goldstone boson (PGB) 
octet in QCD. Spontaneous breakdown of the approximate SU(3) chiral symmetry of QCD by the vacuum, 

SU(3) L x SU(3) R - SU(3) V , (1) 
leads to an octet of PGBs. They are conveniently collected into an SU(3) matrix, 

£ = e^ 21 */^ (2) 

$ = i,j = u 7 d 7 s (3) 

Tr($)=0, (4) 
where / is the leading-order meson decay constant. Under chiral symmetry transformations £ transforms as: 

£ -> (5) 



1 It is important for this argument that the rj' can be integrated out of the PQ theory, which is shown in Ref. 2] . 
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LeSU(3) L , ReSU(3) R (6) 

In the meson sector, xPT is an expansion in powers of e 2 ~ "rn 2 PGB / {AiTf) 2 ~ Pp GS /(47r/) 2 . The lowest-order 
Lagrangian is of 0(e 2 ): 



C 2 = ^Tr(^S^Et) - ^Tr(x£+ + E X ) 



(7) 



Here x is proportional to the quark mass matrix, 

Xij = k^Xi , (8) 
and is defined such that the mass of a PGB containing quarks i and j is 

2 Xi Xj fr.\ 

m l3 = ^ (9) 

at lowest order. Similarly, the 0(e 4 ), or Gasser-Leutwyler, Lagrangian is 2 

£ 4 = ~Li[Tr(d fl T,d fl 'S^)] 2 — L2 r Tr(d fl Y,d l/ l^)Tr(d l j,'Ed l ,Y,^) — L 3 Tr(9 M E<9 M E^(9,,E(9,,Et) 

+i 4 Tr(9 Al Ea M Et)Tr( X Et + E X ) + L 5 Tr(d^d^( X ^ + £ X )) - i 6 [Tr( X E t + S X )] 2 
-L 7 [Tr(xS t - Ex)] 2 - L 8 Tr(xS t xE t + Sx^x) + contact terms . (10) 

It is important to recall that SU(3) group relations and 0{e 2 ) equations of motion were used to reduce the operators 
in £4 to a minimal, linearly-independent basis. 

Partial quenching introduces unphysical, bosonic ghost quarks partners of equal mass for the valence quarks Q: 

E = e^/^ (11) 

$=($ y ), i,j = A,B,..., 1,2,..., A,B,... (12) 

Str($) = (13) 

"Str" indicates supertrace, the graded analog of the trace. The transformation properties of E are the same as in 
Eq. Q , except that the chiral symmetry group becomes a graded group with both commutation and anti-commutation 
relations: 

SU(N sea ) -> SU(N valence 

H~ ^sea I -^valence 

At leading order, the PQ chiral Lagrangian takes the same form as for QCD, Eq. Q, except that traces are replaced 
with supertraces. There are no additional operators invariant under the new PQ chiral symmetry group. It was 
asserted in Refs. 0, 0, that the same is true for the NLO Lagrangian, £4. This is in fact incorrect. The number 
of four-derivative operators in Eq. IjlOl) was reduced from four to three using SU(3) group relations. This reduction 
does not occur in the PQ theory, as we now explain. 

We can use group theory to determine the number of linearly independent four-derivative operators in £4 in a 
general SU(N\M) PQ theory. 3 It is useful to build such terms out the right-handed Lie derivative 4 E^c^E. Under 
the chiral group this transforms as 

(£t0 M £) — > R fit , (15) 



2 We do not consider the Wess-Zumino-Witten Lagrangian in this paper, since it does not contribute to the processes we study until 
higher order than we consider. 

3 In standard discussions for SU(N), the analysis is usually based on the Cayley-Hamilton theorem, from which one can determine 
relations between products of traces of matrices. This approach does not, however, generalize to graded groups, so we use an alternative 
method which does generalize. 

4 All terms involving derivatives in £4, Eq. I1UI . can be written in terms of Lie derivatives by appropriately inserting factors of S'E = 1. 
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FIG. 2: Symmetric product of two adjoints in SU(N\M). The undotted boxes are fundamental representations and the dotted 
boxes are anti- fundamental representations. The result also applies to SU(N) for TV > 3 (with the dashed lines removed). 

and thus potentially contains an adjoint and a singlet component. The singlet is absent because the Lie derivative 
has vanishing supertrace. We can write it as 

Ti%H = U IA ($(x))'?, (16) 

where r is an SU(N\M) generator, and the vector Z„ is a real function of $ and its derivatives. The four-derivative 
terms contain four Lie derivatives, with the Lorentz indices necessarily contracted in two pairs, so that the terms have 
the schematic form 

{^■7){h-r){h-r){h-r). (17) 

Here group indices are not yet contracted, and the issue is to determine the number of singlets under the chiral group 
that are contained in this quantity. The product (l^ ■ f)(l^ ■ f) contains all representations arising in the symmetric 
product of two adjoints. Each of these can be combined to make a singlet with the corresponding representation 
coming from the other product with fj, — ► v. We conclude that the number of independent singlets that one can make 
out of right-handed Lie derivatives is given by the number of representations contained in the symmetric product of 
two adjoints. In fact, this gives all the independent four- derivative terms because all invariants that can be built out 
of four left-handed, or two left-handed and two right-handed, Lie derivatives can be rewritten in terms of right-handed 
Lie derivatives alone, using the cyclicity of the supertrace. 

figure |21 shows the symmetric product of two adjoints in SU(N\M). We use the notation for graded Young tableaux 
of Refs. [a, @ ■ There are four representations, all of which are self-conjugate. They can each form a flavor singlet 
with the same representation in the other pair of adjoints, so there are four linearly-independent operators in the PQ 
theory. Without the dashed diagonal lines, Figgis also the symmetric product of two adjoints in SU(N) for N > 3. 
However, in SU(3), this product contains only three representations: 

(8 ® 8)Symm. = 27 © 8 © 1 (18) 

We conclude that there is one additional four-derivative operator in the PQ chiral Lagrangian, just as there is an 
additional operator in unquenched theories with four or more flavors. 

It is straightforward to determine the form of the new operator. In addition to the three in Eq. (|10f> with trace 
replaced by supertrace, we have: 

Str(0 M E&,Et0"E0"Et) (19) 

For comparison with unquenched QCD, it is convenient to use a linear combination of operators that vanish in the 
unquenched SU(3) sector of the PQ theory: 5 ' 6 

PQ = {str^E^E^E^Et) _ Istr(0 / ,E0"£ t ) 2 

-Str^E&Et) . Str^E^Et) + 2Str(<9 M £<^ I £ t <9„E<9 ,y E t )} . (20) 

The operator Opq should be added to the PQ chiral Lagrangian multiplied by an undetermined coefficient, which we 
call Lpq. This new operator is the subtlety overlooked in Refs. 0,0- 



5 The fact that this combination vanishes in SU(3) can be most easily seen using the Cayley-Hamilton theorem, as nicely addressed in 
Ref. Q. 

6 To be precise, Opq vanishes whenever £ has a block-diagonal form with one block, contained entirely within the ungraded SU(N), 
being an SU(3) matrix, and the other being the identity. 
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FIG. 3: Scattering process to which Ova does not contribute. It is physical when A and B are sea quarks, and unphysical 
when they are valence quarks (with masses differing from those of sea quarks). 
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FIG. 4: Unphysical scattering process because it involves more than three quarks. Here A, B, C, and D are all different. 

III. CONSEQUENCES OF THE NEW OPERATOR AT NLO 

We now discuss the consequences of the additional operator. In general, its presence implies that PQxPT has an 
additional NLO coefficient that affects mesonic quantities, and which must be accounted for in lattice fits. Thus it 
complicates, but does not invalidate, the program of using PQ^PT to extract the real xPT coefficients. In particular, 
since the operator has four derivatives it does not contribute to PGB masses and decay constants, nor to the hairpin 
vertex of 0, at tree-level. It only contributes to these quantities at one-loop order, which is NNLO in the chiral 
expansion. We return to these contributions in the next section. 

In this section we consider the quantities to which Opq does contribute at NLO. These are tree-level scattering 
processes, involving at least four PGBs. We focus on the simplest, four-meson, scattering processes. We illustrate 
such processes using quark line diagrams, which show the flow of quarks within the mesons. 7 Opq contains both 
single and double supertraces, so it contributes to both connected and disconnected quark line diagrams. Here we 
discuss the types of charged meson scattering processes to which Opq does and does not contribute. 8 

Because O-pq vanishes in the unquenched SU(3) limit, it cannot contribute at tree-level to scattering processes 
that occur in the unquenched SU (3) sector. In fact, since the operator has no dependence on the quark mass (i.e. no 
factors of x), it cannot distinguish between sea and valence quarks. In other words, its contributions are the same as 
in the chiral limit, in which limit there is a full SU(N va i ence + N sea ) symmetry among quarks. It follows that any 
scattering process involving three or fewer quarks, whether they are valence or sea quarks, is equivalent to an SU (3) 
process, and cannot receive a contribution from Opq. One such process is shown in Fig. [21 Mathematically, the 
contribution of Opq to this process vanishes because relative minus signs and numbers of Wick contractions among 
the quark line diagrams produce cancellations. Note that these cancellations would not occur if either or both of the 
quarks in Fig. |3 were bosonic. 

Opq does contribute to scattering processes that are unique to the PQ theory, and thus unphysical. For example, 
that shown in Fig. 0] involves four quark flavors. In the case of PQ QCD with N sea = 3, these could be four valence 
quarks, three sea quarks and a valence quark, or other possibilities, figure is also unique to the PQ theory, involving 
both valence and ghost quarks. It is clear that Opg contributes to these processes because there is only one quark- line 



7 Quark lines follow the flavor indices of the meson fields within each supertrace. 

8 We do not consider processes involving neutral mesons because double poles in the meson propagators make it unclear how to amputate 
diagrams and thus define scattering amplitudes even at tree level. 
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FIG. 5: Unphysical scattering process because it involves ghost quarks. 



diagram and thus no possibility of cancellations. These examples illustrate the origin of the additional operator: the 
extra particles present in the PQ theory allow the separation of individual quark contractions in a way not possible 
in the unquenched sector. The same explanation holds for physical theories with four or more flavors. 

These unphysical scattering processes can be used, in principle, to determine Lpq, the coefficient of the new operator 
in the PQ chiral Lagrangian. To illustrate this we have calculated some representative scattering amplitudes. First, 
we give the expression for the "physical" process of Fig. [3] 
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(21) 



where the momenta are Euclidean and are defined in the figure, and m A B is the leading order mass as given in Eq. 0. 
Note that, as claimed above, the amplitude does not contain a term proportional to Lpq. This is true even though 
the process may involve only valence quarks with masses differing from those of sea quarks, and thus be unphysical. 

The unphysical processes in Figs. 21 and El are chosen so that they only receive contributions from disconnected 
quark line diagrams. It follows that their scattering amplitudes only receive contributions from double supertrace 
operators, and both have terms proportional to Lpq. The amplitude for Fig. 0] is 



Ma = 



-ji L im ■P3)(P2 - Pa) - jgL 
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16 r 

-j4 L PQ[((Pi -P2){P3 - Pa) + (pi -P3){P2 - Pa) + (Pi -Pa)(P2 ■ 

16 , r 2 / \ 2 / ^ ^ 59 

j^L^mcjjipi ■P3) + m AB {p 2 - Pa) 

— T (sT I ( m AB + m cp) 
fi 6 I / i % sea I 2 ' 

J \sea / 



- JZ L e m AB m CD, 



(22) 



The amplitude for Fig.0is really a special case of Eq. H22J) . obtained by setting xc = Xa and xd = Xb, and including 
an overall sign because Fig. [S] contains one pair of bosonic quarks, and therefore one supertrace over bosonic indices. 
The kinematic factors are identical because the quark line diagrams have the same structure. 

These amplitudes provide a method, in principle, of determining Lpq from PQ simulations. The idea is to calculate 
the scattering amplitudes, as a function of the valence masses and momenta, and thus determine Lis horn M.3 
(£4,5,6,8 having been obtained from fits to the PGB masses at NLO Q, 0), and then determine Lpq from VW4. 
Other amplitudes involving different external quarks can also be used. In practice, it is notpossible to determine the 
scattering amplitude itself, because the standard method for doing so relies on unitarity y|, and this is violated at 
the one-loop level in these amplitudes Instead, what must be done is to use PQxPT to predict the form of 

Euclidean finite volume correlation functions, which will be given in terms of Lpq and other low energy constants, 
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FIG. 6: PGB mass renormalization from Oj>q (represented by the solid circle). V and S refer to valence and sea quarks 
respectively. 
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FIG. 7: PGB decay constant renormalization from O-pQ- The filled square represents the axial current, while the filled circle 
represents interactions from the chiral Lagrangian. The label NLO indicates where Ovq is contributing. 

and fit this to lattice results for these correlation functions. 9 In this way the unphysical nature of the PQ theory is 
accounted for. 

IV. CONSEQUENCES OF THE NEW OPERATOR AT NNLO 

As noted above, O-pQ contributes to PGB masses and decay constants first at NNLO in PQxPT. We have calculated 
these contributions as a first step in extending the application of PQ QCD to this order. To simplify calculations we 
consider only two valence quarks, A and B, and N sea quarks of equal mass. We also only consider the properties 
of charged mesons, ttab- figures [S] and {7\ show the diagrams contributing to the mass and decay constant of charged 
mesons. All of the meson diagrams receive contributions from the classes of quark line diagrams shown in figure |S] 
Quark line diagrams show the effect of partial quenching — valence quarks appear only as external states because ghost 
quarks cancel them in loops. The series of sea quark hairpin diagrams removes the flavor singlet propagator, enforcing 
the fact that graded generators are traceless. 

The corrections to the renormalized meson mass and decay constant from the new PQ operator (using dimensional 
regularization and the MS scheme of Ref. are: 
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The simplest example of this is the hairpin correlator discussed in Ref. which is predicted to have a double pole, and thus is manifestly 
unphysical. Nevertheless, its coefficient, if it can be determined, gives information about physical low energy constants L5 and L7. We 
note that in the quenched theory, where one also expects a double pole, its coefficient has been successfully determined from fits. 
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FIG. 8: Quark line diagrams contributing to mass and decay constant renormalization. 
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FIG. 9: NNLO PGB mass renormalization from operators in C% 
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(24) 



where /i is the renormalization scale. It can be seen that the corrections in Eas. 1231 and 1241 are NNLO by the fact that 
they have the form x 2 ln(x). A check of these results is that if A and B are sea quarks, and if N sea = 3, the mesons 
in Figs. El and live in the unquenched SU(3) sector and cannot get contributions from Opq. Equations l)23|l and 
(|24|l do, in fact, vanish in the unquenched limit: \A — Xb — Xs, N sea = 3. 10 

We stress that these results are only part of the NNLO charged meson mass and decay constant corrections. The 
full correction receives 2-loop contributions from C2, 1-loop contributions from terms in £4 other than O-pq, and 
tree-level contributions from Cq. We next discuss the contributions from £q. 



V. ANALYTIC NNLO MASS AND DECAY CONSTANT CORRECTIONS 

Fits of present PQ lattice data require NNLO terms Full non-analytic NNLO calculations in PQxPT 

are not available, so, as an intermediate step, we have determined the form of the analytic NNLO mass and decay 
constant corrections for charged mesons. These arise from the tree-level diagram in Fig. 

This diagram gets contributions from operators in £ 6 with zero and two derivatives. Analytic, NNLO decay constant 
corrections come from a similar diagram, but with an axial current insertion, and only come from two-derivative 
operators. As in the case of Opq, we need to be sure to include all operators in the PQ theory, some of which may 
be linearly-dependent in unquenched SU{3). 



The same holds true for N aea = 2, where O-pQ also vanishes. 
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First we determine the number of linearly-independent, derivative-free operators in £§. As always, these must be 
invariant under chiral symmetry transformations, hermitian, and invariant under the discrete symmetries respected 
by QCD: C, P, and T. They must have three x matrices to be 0(e 6 ). It is easy to construct terms with the correct 
symmetry properties using the following two building blocks with simple transformation properties: 

E+X - {Rx}tf)(L X I&) = i?(E t x)i? t (25) 

We find seven allowed forms: 

1. (str( X t S )) 3 +(str(£t x )) 3 

2. Str^Ex^xtE) + Str(Et xS t xS t x ) 

3. Str(x t Ex t S)Str(x t E) + Str(StxE t x)Str(S t x) 

4. Str(x t E X t E)Str(Et x ) + Str(EtxEt x )Str(x t E) 

5. (str( X tE)) 2 Str(Et x )+ (str(Et x )) 2 Str(x t E) 

6. Str(x t EEt x )Str(x t E) + Str(x t EEt x )Str(Etx) 

7. Str(x t EEt x Etx) + Str^xx^E) 

(The last two can be simplified using EtE = 1, but we choose not to do so to better show their form.) Our task is to 
reduce them to a linearly independent set using group-specific relations. 

For SU (3) we can use the Cayley-Hamilton theorem, which states that any matrix satisfies its own characteristic 
equation, and gives rise to numerous relations between traces of matrices. The relevant relation here is that the 
determinant of a matrix can be expressed as a polynomial in traces of the matrix of order the size of the matrix. For 
3x3 matrices, this polynomial is 

det(M) = -Tr(M 3 ) - -Tr(M 2 )Tr(M) + -(Tr(M)) (26) 

Since det(x^S) = det(x)* = const., the first three operators in the list are linearly dependent. Thus there are only 
six independent derivative-free 0(e 6 ) operators in SU(2>). 

There is no analogous relationship among operators in the PQ theory, because the super-determinant is not a 
finite polynomial 0- The same holds true for SU(N), N > 4, since the Cayley-Hamilton relations involve quartic or 
higher-order polynomials, and do not imply any linear dependence among third order polynomials. We conclude there 
are seven derivative-free operators in PQxPT, and also in unqucnched SU (N) theories with N > 4. One of the seven 
PQ operators can be chosen to be unphysical, namely the linear combination given in Eq. Ij26|l with trace replaced 
with supertrace. While our PQ result is new, we note that Ref. ^| enumerates all of the operators in the C(e 6 ) chiral 
Lagrangian for unquenched SU(N). Ref. [lj lists seven linearly-independent operators in SU(N), N > 4, and six in 
SU(3), so our results our consistent. 

For terms involving two derivatives, we can use group theory to determine the number of independent operators, 
generalizing the method used earlier in Sec. |nj To be 0(e 6 ) such operators must also contain two factors of \. 
We construct these operators using right-handed Lie derivatives and the matrices defined in Eq. (|25[) . all of which 
transform non-trivially only under the right-handed chiral group. We recall from Sec. ^ that the Lie derivative 
transforms like an adjoint, so an operator with two Lie derivatives (which are necessarily contracted together by 
Lorentz symmetry) must be symmetric in this pair of adjoints. The symmetric product of two adjoints is shown in 
Fig. [5] E^x a- n d X ^ transform as bi- fundamentals under the chiral symmetry group, but are not traceless, and so 
contain both adjoint and singlet parts. Therefore any operator that contains two of these matrices also comes from 
the product of an adjoint plus singlet times an adjoint plus singlet, as shown in Fig. 1101 If the operator contains two 
of the same matrices (either E^x or ^'E), it will only come from the symmetric part of this product, shown in Fig. II II 
Putting this together, we must construct singlets out of the product in Fig. [3 which comes from the Lie derivatives, 
and that in either Fig. I1UI or II II which come from the terms involving factors of x- 

In this way, we find that eight linearly-independent operators can be made out of two Lie derivatives, one E^x, and 
one x^- They correspond to the eight flavor singlets in the direct product of the representations in Figs. l2l and 1101 
One choice of basis is the following, 

1. Str(5 AI Et5 AI Ex t x) 
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FIG. 10: The product of two bi-fundamentals in SU(N\M). The result also applies for SU(N) for N > 3 (with the dashed 
lines removed). 



(H3 + 1) >< (HQ + 1) = ( HELZH + 

^ ' ^ ' Symm. ^ 

FIG. 11: Symmetric product of two bi-fundamentals in SU(N\M). The result also applies to SU(N) for TV > 3 (with the 
dashed lines removed). 



2. Str(5 M St xx ta /i E) 

3. Str(^£t xS t^ Sx t S ) 

4. Str^Et^strO^x) 

5. Str(^St x )Str( X t^S) 

6. Str(5 AI Eta AiS] )Str(x t S)Str(Et x ) 

7. Str(5 AI Eta AiS]x t S )str(Et x ) 

8. Str(5 AI Eta AiES] t x )str( X tE) 

where the derivatives act only on the quantity immediately to their right, and we have used the anti-Hermiticity of 
the Lie derivatives to simplify some of the operators. Although group theory gives operators with the right chiral 
transformation properties, we must impose the other symmetries by hand. Operators 1 and 2 transform into each 
other under parity, so we must include their sum, with a single coefficient, in the Lagrangian. Operators 7 and 8 
are hermitian conjugates, so we must include their sum as a single operator as well. We conclude that there are six 
operators of the above type in the PQ chiral Lagrangian. 

Similarly, six operators can be made out of two Lie derivatives and two £^x s or X^ s because there are two fewer 
adjoints in the symmetric product of Fig. 1111 

1. Str^St^Sx+Sx^ + d^d^x^X)) 

2. Str(^£t x d M £t x + x t^Ext^£) 

3. Str(5 M St^S)Str(x t Sx t S + ^x^x) 

4. Str(x t 9 AI S)Str(x t ^ l S) + Str(^Et x )Str(^St x ) 

5. Str(a Al Eta Al I]x t S)Str(x t S) + Str^Et^SSt^Str^^) 

6. Str(a AI Eta AiS] )str(x t S)Str(x t S) + Str(^£t^£)Str(£t x )Str(Et x ) 

This completes the counting for the PQ theory: there are twelve linearly independent two-derivative operators in 
Cq. The same result applies to SU(N) for N > 3. For SU(3), however, one operator can be eliminated from each 
of the two previous lists in because there is one fewer representation in the product of adjoints for SU (3) than for 
SU(N\M). We can use Cayley-Hamilton relations analogous to Eq. (|26H to obtain the precise relationship between 
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the operators. 11 Thus there are only ten two-derivative operators in SU(3), and correspondingly two unphysical 
operators in the PQ theory. Our results for SU(N) are consistent with those of Ref. [l5j . 

We stress that the counting of operators in all SU(N\M) PQ theories is valid for any N > M > 2. The super- 
determinant is not a finite polynomial in any such theory, and there are no Cay ley-Hamilton relations. From the 
perspective of Young tableaux, this result follows because the theories do not have an antisymmetric tensor to 
contract indices, and in particular there is no limit to the number of boxes in the column of a Young tableau, as there 
is in SU{N). 

We now return to the NNLO corrections to the PGB mass and decay constant. Although there are contributions 
from all of the operators in Ce enumerated above, the form of the tree- level correction to the mass of ttab is simple: 12 



Sm AB 

. TO 2 , D / NNLO, 



"AB 



Analytic 



ai^Us) + a z tT (xs) 2 + +a 3 te(xs)(xA + XbJ + a 4 (xA+ Xb^J + «5 (xa - Xb) • (27) 



Here xs is the mass matrix of the sea quarks. In fact, Eq. I|27|) is the most general quadratic polynomial separately 
symmetric in the valence and sea quark masses. The cti are linear combinations of Cq coefficients, and we have checked 
that they are independent, so that no relations are predicted among the cti. Corrections to Jab have the same form 
with different coefficients. Fits in Ref. 0] to PQ lattice data with \mg < m sea < |to_5 are consistent with this 
NNLO formula. 

It is interesting to consider the contributions to the on from the three unphysical operators. We find that these 
operators do contribute to the cti (in fact, to all five of them in the basis we use), i.e. that the unphysical operators 
affect PGB properties at tree level. 13 This is different from Opq, which does not affect PGB properties until one-loop 
order. It implies that the on are linear combinations of the coefficients of physical and unphysical operators. Clearly, 
since the number of low energy constants entering at NNLO (nineteen in the zero and two-derivative sector alone) 
exceeds the number of constraints from PGB masses and decay constants (ten in the PQ theory 14 ), one must use 
other quantities in order to completely determine all the constants. 

Just as at NLO 3], using the PQ theory can simplify this determination. Consider the case of degenerate sea 
quarks. The unqucnched sector alone constrains only two combinations of constants (one each from the single term 
proportional to x\ = Xb = Xs m ^ m AB an( i ^Iab), although these combinations necessarily involve only physical 
constants. Consideration of the PQ theory adds three unphysical operators, but this "cost" is outweighed by the 
benefit of the six additional constraints from Eq. (|27J) . 15 In addition, using the PQ theory one can add more values 
of \a and xb to the fit at relatively small computational expense. 

If one uses non-degenerate sea quarks then only the second benefit of partial quenching applies. This is because, 
even in the unquenched sector, one can, in principle, determine four of the five on for both PGB masses and decay- 
constants, using a number of choices for the sea quark masses. 16 Moving to the PQ theory gives two more constraints 
but at the cost of three unphysical operators. Nevertheless, we suspect that this cost will be easily outweighed by the 
additional data points that one can obtain in the PQ theory. 



VI. CONCLUSION 



In this paper we have discussed a complication to the program of using simulations of unphysical PQ theories to 
obtain physical low energy constants. We find that there are operators in PQxPT which vanish when restricted to 
the physical, unquenched sector of the theory. There is one such unphysical operator in the NLO Lagrangian £4, 
and three in the zero and two derivative part of £q. Generically, we expect the number of unphysical operators 
to increase with the order of the calculation. Their presence is related to the fact that correlation functions in the 



11 See Ref. for details. 

12 This result has also been obtained independently by Ref. Il4l . 

13 Of course, by definition, the total contribution of unphysical operators to 5m Ag must vanish when both \A an d \B are equal to sea 
quark masses, and we have checked this explicitly for the combination appearing in Eq. with N sea = 3. 

14 There are five each from Eq. 1271 and the analagous expression for S/ab ■ 

15 There are three each from 6m ? AB and S/ab > rather than four, because the first two terms collapse into a single term for degenerate sea 
quarks. 

16 With three sea quarks there are four independent quadratic combinations of masses that are symmetric under A <-> B: {xa ±Xs) 2 > Xq 
and (\A + Xb)XC- 
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unqucnched sector involve linear combinations of quark contractions, while in the larger PQ theory one can determine 
each contraction separately. 

Dealing with this complication is, in principle, straightforward. The contributions from unphysical operators can 
be disentangled from those of physical operators as long as one considers enough quantities in the fits. The extent 
of this complication will depend on the example being considered. For the quantities we have studied in detail, the 
PGB masses and decay constants, the new C(p 4 ) operator Opq does not contribute at tree level to PGB masses and 
decay constants, and so one does not need to include it at all until one works at NNLO. 

As a case study in how to deal with the complication of unphysical operators, we have outlined a strategy for how 
one could, in principle, determine and subtract the NNLO contributions of Opq to PGB properties. Our idea is to 
determine its coefficient using unphysical scattering processes, to which Opq contributes at the first non-vanishing 
order, which is NLO in xPT. Then, when doing a NNLO fit to the PGB properties, one can subtract the contribution 
from Opq, the form of which we have calculated. Of course, any such fit would require all other NNLO terms, 
including two-loop contributions with vertices from £2, and these are not yet available for the PQ theory. 

We have also considered the more difficult question of how to separate the contributions from the unphysical 
operators in C§. The difference here is that these operators contribute to PGB masses and decay constants at tree 
level, so their coefficients appear in combination with those of the physical NNLO operators. One must therefore rely 
on the general strategy of calculating enough quantities to constrain all the coefficients, both physical and unphysical. 

Clearly, extending the determination of low energy constants to NNLO in the meson sector is a significant under- 
taking, requiring extensive x?T calculations beyond those presently available, and consideration of many quantities in 
addition to PGB masses and decay constants. If such a program is undertaken, PQ simulations can play an important 
role. While they introduce a few additional unknown constants, they provide both additional constraints on low 
energy constants from the extra functions of valence masses that are available, and, most importantly, the possibility 
of additional, relatively cheap, data points to include in the fits. 

Finally, we take this opportunity to reiterate the purpose of simulating PQ QCD. The intent is not to study it as 
an approximation to, or a model of, QCD, for it is an unphysical theory, but rather to use it as a tool to more easily 
extract the low energy constants of QCD. Once these are determined, they can be used to calculate PGB processes 
in QCD, including those that are not easily accessible to lattice calculations. The unphysical nature of PQ QCD is 
manifested in several ways, for example by the presence of double poles in neutral propagators, and by the inability to 
define decay and scattering amplitudes pJJ |. We have discussed here another unphysical feature, namely the presence 
of additional operators in PQxPT which do not contribute in the unquenched sector. 17 While all these features do 
complicate the use of PQ QCD, we have argued that these complications can be overcome and that PQ simulations 
remain an important element of the lattice practitioner's toolkit. 
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